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ABASTRACT 

The  objeetive  of  this  researeh  is  to  develop  a 
modeling,  simulation,  and  design  tool  for  nano- 
eomposites  in  order  to  prediet,  with  aeeepted  fidelity, 
nano-eomposite  behaviors;  and  to  design  optimum  nano- 
eomposites  for  speeifie  Army  applieations  sueh  as 
ballistie  proteetion.  Based  on  a  unit  eell  model  ereated 
for  nanoelay-epoxy  eomposites,  the  effeet  of 
nanopartiele  distribution  on  the  maximum  stress 
developed  in  epoxy  resin  was  investigated  using  the 
Meshfree  Partiele  Method  based  simulation  tool 
developed  at  MKP.  The  ensemble  average  of 
meehanieal  property  for  nanoelay-epoxy  eomposite  was 
also  studied.  The  inerease  in  the  meehanieal  properties 
for  nanoelay-reinforeed  epoxy  was  evideneed  from  the 
simulations.  It  is  found  that  statistieal  analysis  is  emeial 
to  prediet  the  performanee  of  nanopartiele  reinforeed 
eomposite.  Several  armor  samples  reinforeed  with 
nanoelay-epoxy  eomposites  were  fabrieated  and  ballistie 
tests  were  eondueted  on  the  nanoelay-reinforeed 
eomposite  armors.  It  was  found  that  higher  meehanieal 
properties,  sueh  as  modulus  and  strength,  will  help  to 
improve  the  performanees  of  eomposite  armor  system. 
A  numerieal  simulation  was  initiated  to  investigate  the 
proteetion  differenee  between  armor  with  nanoelay 
reinforeed  epoxy  and  armor  with  pure  epoxy  resin. 
Efforts  have  been  made  to  reproduee  similar  ballistie 
proteetion  behavior  using  the  simulation  tool  in 
eomparing  with  the  aetual  test  results. 

1.  INTRODUCTION 

Inspired  by  the  exeellent  results  obtained  by  a  group 
of  researehers  in  Japan  by  the  embedment  of  nanoelay 
into  nylon  [1],  extensive  work  has  been  earried  out  over 


the  past  deeades  on  the  ineorporation  of  nanoelay  into 
polymers  for  the  purpose  of  enhaneing  the  polymer 
properties.  Normally  the  properties  under  investigation 
inelude  tensile  modulus,  tensile  strength,  glass  transition 
temperature,  resistanee  against  the  absorption  of 
moisture,  resistanee  against  flammability,  and  fraeture 
toughness.  Ahmadi  et  al.  [2]  used  X-ray  diffraetion  and 
transmission  eleetron  mieroseopy  to  investigate  the 
morphology  of  elayed-reinforeed  ethylene  propylene 
diene  terpolymer  (EPDM)  nanoeomposites  with  organo- 
montmorillonite  (OMMT).  It  was  found  that  the  tensile 
modulus  ean  be  inereased  by  60%  and  tensile  strength 
ean  be  inereased  by  175%.  Avila  et  al.  [3]  investigated 
the  penetration  meehanism  on  polymer-nanoelay-fiber 
glass  nanoeomposite  using  low-veloeity  impaet  tests 
(ASTM  D5628-01).  The  sample  was  made  of  S2-glass 
plain-weave  woven  fabrie  and  nanoelay-epoxy.  It  was 
shown  that  the  front  side  delamination  area  was  redueed 
22%  with  1  weight  %  nanoelay  addition.  The  baek  area 
de lamination  reduetion  ranged  from  10-21%  with  respeet 
to  different  impaet  energy.  It  was  also  found  that  slow 
eure  in  epoxy-nanoelay  mixture  leads  to  mueh  better 
nanoeomposites  without  bubble  formation.  Wetzel  et  al. 
[4]  used  a  Charpy  impaet  strength  test  to  investigate  the 
fraeture  toughness  of  nanoeomposites.  The  test  sample 
is  eomposed  of  epoxy  resin,  titanium  dioxide  (Ti02) 
nano-partieles  and  ealeium  silieate  (CaSiOs)  miero- 
partieles.  The  improvement  in  stiffness,  impaet  strength 
and  wear  resistanee  were  witnessed  at  low  nano-partiele 
filler  eontents.  Optimal  Ti02  nano-partiele  eontent  is 
found  at  4  volume  %.  Additional  3  volume  %  CaSiOs 
miero-partieles  into  epoxy/Ti02  (4  volume  %)  ean 
further  improve  the  impaet  toughness  of  the 
nanoeomposite.  Yao  et  al.  [5]  earried  out  dynamie  tests 
on  a  nanoelay-epoxy  beam  with  a  manufaetured  initial 
eraek.  The  beam  is  made  of  3  weight  %  Montmorillonite 
(MMT)  silieate  elay  partieles  and  epoxy  resin.  It  was 
found  that  fraeture  toughness  improved  with  the  addition 
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of  clay  particles.  The  toughness  improvement  can  be 
attributed  to  the  increasing  interfaces  from  particle 
reinforcement  at  the  nanoscale.  Zheng  et  al.  [6]  carried 
out  impact  tests  (ASTM-D256)  and  tensile  tests  (ASTM- 
D638)  for  nanoclay  (Si02)  -  epoxy  composites.  For  a 
glass  fiber  nanocomposite,  additional  three-point 
bending  tests  (ASTM-D790),  compression  tests  (ASTM- 
D6641)  and  shear  tests  (ASTM-D2344)  were  also 
carried  out.  Two  kinds  of  nanocomposite  samples  were 
made  for  the  mechanical  tests:  nano-particle  (Si02)- 
epoxy  composite  and  glass-fiber-epoxy  nanocomposite. 
For  nano-epoxy  polymer,  it  was  found  the  tensile 
strength  increased  115%,  tensile  modulus  increased  13% 
and  impact  strength  increased  60%  at  the  optimal  3 
weight  %  of  nano-particle  addition.  From  the  impact 
test,  it  was  seen  that  small  curvatures  occurred  at  the 
surface  of  the  nanocomposite,  indicating  that  massive 
new  surfaces  were  created  under  impact  loading.  As  a 
result,  much  more  impact  energy  could  be  dissipated 
through  these  new  surfaces  in  the  Si02-epoxy 
nanocomposite.  For  the  nano-glass  fiber  composite,  the 
mechanical  property  improvements  were  found  as: 
bending  strength  69.4%,  tension  modulus  21%,  tensile 
strength  23%  at  5  weight  %  nano-particle  addition. 
Flowever,  shear  and  compression  strength  enhancement 
were  limited.  Fan  et  al.  [7]  described  various  different 
nano-polymer  applications,  especially  for  nylon  6 
nanocomposite  with  2  weight  %  nanomer  and  for 
Polyolefin  nanocomposites  (FIPP).  For  nylon  6 
nanocomposite,  the  property  enhancement  reached  its 
optimal  value  at  2  weight  %  nano-particle  addition  and 
the  improvement  was  reported  as:  flexural  modulus  35%, 
tensile  modulus  28%.  For  FIPP  with  6  weight  %  nano¬ 
particle,  the  property  enhancement  was  even  greater: 
98%  in  tensile  modulus  and  78%  in  flexural  modulus. 

The  addition  of  nano-particle  is  also  beneficial  to 
composite  chemical  and  electrical  properties.  Bharadwaj 
[8]  investigated  the  relative  permeability  of  layered 
nanoclay  samples  based  on  a  modified  tortuosity  model. 
The  effects  of  sheet  length,  concentration,  orientation, 
and  degree  of  delamination  on  the  relative  permeability 
were  explored.  The  layered  silicate  clay  was  1  nm  in 
thickness.  The  length  of  the  sheet  ranged  from  30  - 
2000  nm.  It  was  found  that  the  gas  barrier  property  was 
enhanced  with  less  than  5  weight  %  nanoclay  addition. 
Pan  et  al.  [9]  investigated  electricity  conductivity  for  one 
kind  of  nanocomposite  (nylon  6  with  0.5 -2. 5  weight  % 
nano  graphite  flake).  The  nano  graphite  was  shown  to 
be  a  good  conducting  filler  for  polymers.  Singh  and 
Balazs  [10]  discovered  that  grafting  short-chain 
surfactants  to  the  clay  sheets  could  increase  the  gallery 
spacing  and  promote  polymer  permeation.  Self 
consistent  field  (SCF)  theory  was  used  to  investigate  the 
effect  of  varying  the  properties  of  both  the  surfactants 
and  polymers  on  the  miscibility  of  the  polymer/clay 
mixture.  It  was  found  that  the  blending  of  polymers  and 


clay  under  optimal  conditions  could  yield 
nanocomposites  with  a  great  tensile  strength,  heat 
resistance  and  gas  permeability  compared  to  the  pure 
polymer  matrix.  Furthermore,  only  1-10  weight  %  nano 
additives  were  needed  for  performance  enhancement. 
For  detailed  discussion  of  chemical  and  electrical 
properties  of  nanocomposites,  the  review  paper  by 
Flussain  et  al.  [11]  is  suggested. 

Numerical  simulation  for  nanocomposites  has  also 
been  a  research  focus  in  recent  years.  Flackett  et  al.  [12] 
used  Monte  Carlo  and  molecular  dynamics  (MD) 
computer  simulations  to  explore  the  atomic  scale 
structure  and  dynamics  of  intercalated  PEO-layered 
silicate  nanocomposites.  Srivastava  et  al.  [13]  also 
studied  the  computational  nanotechnology  with 
molecular  dynamics  for  fullerenes  and  carbon  nanotube 
(CNT)  based  molecular  materials.  It  has  been 
emphasized  that  quantum  MD  is  needed  for  the 
simulation  in  picosecond  (10"'^  second)  scale,  while 
classic  MD  is  commonly  used  to  the  situation  where  the 
motion  of  atoms  and  molecules  is  treated  using  finite 
difference  equations  of  Newtonian  mechanics.  Tight- 
binding  molecular  dynamics  (TBMD)  was  proposed  as  a 
bridge  between  quantum  MD  and  classic  MD  [13]. 

Instead  of  molecular-dynamics-based  simulation 
method,  Li  and  Chou  [14]  modeled  the  deformation  of 
CNT  using  a  structural  mechanics  approach.  Potential 
energy  was  used  to  model  the  bonding  between 
molecules.  The  tensile  modulus  and  shear  modulus  of 
single-walled  carbon  nanotube  were  examined.  Zhao 
and  Floa  [15]  also  examined  tensile  modulus,  tensile 
strength  and  fracture  toughness  of  a  nanoclay  reinforced 
epoxy  cell  using  FEM  simulation.  The  amount  of 
nanoclay  addition  modeled  was  5  volume  %.  It  was 
shown  that  nano-particle  volume  percentage  is  the  main 
factor  for  increased  elastic  modulus.  Stress 
concentration  in  nanocomposite  depends  highly  on  the 
nano-particle  dispersion.  The  more  dispersive  the  nano¬ 
particle  distribution,  the  less  stress  concentration  exists 
in  the  composite. 

The  simulation  techniques  can  be  summarized  as 
follows.  In  principle,  any  problem  associated  with 
molecular/atomic  motions  can  be  simulated  by  molecular 
dynamics,  tight-binding  molecular  dynamics  and  density 
function  theory  (DFT).  Flowever,  due  to  their  huge 
computational  tasks,  practical  applications  of  these 
atomistic  modeling  techniques  are  limited  to  systems 
containing  a  small  number  of  molecules  or  atoms  and  are 
usually  confined  to  studies  of  relatively  short-lived 
phenomena,  from  picoseconds  to  nanoseconds. 

The  continuum  mechanics  modeling  approach  can 
be  used  to  analyze  the  static  or  dynamic  mechanical 
properties  of  composites  with  nanoparticles.  Flowever, 
these  models  neglect  the  detailed  characteristics  of  nano¬ 
particle,  and  are  unable  to  account  for  forces  acting  on 
the  individual  atoms.  Therefore,  there  is  a  need  for 
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developing  a  modeling  teehnique  that  analyzes  the 
meehanieal  response  of  nanoeomposites  at  small  seale 
but  is  not  overburdened  by  the  time  seales.  Sueh  a 
modeling  approaeh  would  benefit  novel  nanodeviee 
design  and  multi-seale  simulations  of  nanosystems.  We 
will  extend  the  theory  of  elassie  stmetural  meehanies 
into  the  modeling  of  nano-partiele  ineorporated  polymer. 
The  bonding  between  nano-partieles  will  be 
approximated  as  potential  funetions  of  moleeular 
meehanies.  The  nanoeomposite  property  will  be 
simulated  in  a  miero-seale  unit  eell.  The  obtained 
meehanieal  property  ean  then  be  used  in  maero-seale  for 
design  optimization. 

2.  DEVELOPMENT  OF  SIMULATION  MODEL 

2.1  Nanoclay  Simulation  Model 

A  multiseale  deeomposition  [16]  was  utilized  to 
eouple  simulations  at  different  seales.  The  fundamental 
idea  is  to  deeompose  the  total  displaeement  field  u(x) 
into  eoarse  and  fine  seales  [17] 

u(x)  =  u(x)  +  u'(x) ,  (1) 

where  the  eoarse  seale  U  is  that  part  of  the  solution 
whieh  ean  be  represented  by  a  set  of  basis  funetions,  i.e. 
finite  element  shape  funetions.  The  fine  seale  is 
defined  as  the  part  of  the  solution  whose  projeetion  onto 
the  eoarse  seale  is  zero.  An  example  of  this 
deeomposition  is  to  model  the  epoxy  region  using  a 
eoarse  seale  and  model  the  nanoelay  region  using  a  fine 
seale: 

u  =  Nd  +  Qq.  (2) 

Here,  N  is  the  shape  funetion  set  for  the  eoarse  seale 
domain  and  d  is  the  FEM  nodal  displaeement  veetor.  In 
the  fine  seale  region,  q  is  the  “exaet”  solution  obtained 
from  any  atomie  or  moleeular-level  simulation  tool,  i.e. 
moleeular  dynamies.  Q  is  the  eomplimentary  projeetion 
operator  defined  in  referenee  [17]. 

With  the  definition  in  Eq.  (2),  a  Lagrangian  L  ean  be 
defined  to  be  the  differenee  between  the  kinetie  energy 
and  the  potential  energy, 

I(u,u)  =  ^U^M>-[/(u),  (3) 

where  M^  is  the  mass  matrix  obtained  from  FEM  and 
MD,  U  is  the  interatomie  potential  energy,  and  veloeity 

u  =  Nd  +  Qq.  (4) 

Substituting  Eqs.  (2)  and  (4)  into  the  Lagrangian  (3), 
results  in 

I(d,d,q,q)  =  |d"Md  +  lqXq  (5) 

-f7(d,q) 


where  the  eoarse  seale  mass  matrix  M  is  given  by 
M  =  N^M^N  ,  and  the  fine  seale  mass  matrix  M^  is 


m^=q"m,q 

The  multiseale  equations  of  motion  are  obtained 
from  the  Lagrangian  [17]  and  have  the  following  form: 

Md  =  N^f  ,  (6) 

M^q  =  Q^f,  (7) 

in  whieh  the  interatomie  foree  f  ean  be  expressed  as 

5f7(u) 


f  = 


0U 


(8) 


A  Lennard- Jones  (LJ)  potential  ean  be  used  to 
define  the  interatomie  potential  for  nano  materials: 


f7(x.,x.)  =  f7(r)  =  4^ 


'  <j  ' 


\r  J 


f 

'  <j  ' 


\r  J 


where  r  is  the  interatomie  distanee,  (T  is  the  eollision 
diameter,  and  S  is  the  bonding  energy. 


2.2  Meshfree  Particle  Method 


A  Meshfree  Partiele  Method  [17,  18]  is  applied  to 
the  simulation  of  nanoeomposites.  The  nanoeomposite 
field  is  diseretized  into  partieles  with  desired 
volume/mass.  Compared  to  eonventional  grid-based 
methods,  e.g.  the  finite  element  approaeh  (FEM),  the 
meshfree  method  is  a  more  powerful  way  to  simulate 
dynamie  problems  sueh  as  high-veloeity  impaet.  These 
problems  often  involve  numerieal  diffieulties  sueh  as 
large  deformation,  material  inhomogeneity,  moving 
material  interfaee  and  deformable  boundary,  whieh  make 
grid-based  methods  ehallenging  in  numerieal 
simulations.  Furthermore,  the  nanoeomposite  domain  is 
eonsidered  as  a  set  of  diserete  physieal  partieles  rather 
than  a  eontinuum.  Consequently,  the  meshfree  method 
is  a  natural  ehoiee  in  simulations. 

The  interatomie  foree  on  partiele  i  is  denoted  as  f„ 
whieh  is  the  summation  of  interaetions  from  neighboring 
partieles, 

V, 

f  =  yf , 

y=i 

(10) 

where  fy  is  the  internal  foree  from  partiele  j  to  i,  Nj  is  the 
total  number  of  partieles  in  the  support  domain  of 
partiele  i. 

The  definition  of  support  domain  [18]  of  partiele  i  is 
illustrated  in  Fig.  1.  The  dimension  of  the  support 
domain  for  partiele  i  is  denoted  as  /?,.  A  partiele  j  is 
eonsidered  within  the  support  domain  if  the  distanee 
between  partiele  j  and  i,  Dy,  is  less  than  the  dimension  of 
the  domain,  hi. 
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The  interaction  force  between  particle  i  and  j,  fy, 
depends  on  the  interatomic  potential  defined  in  Eq.  (8). 
For  the  L-J  potential  defined  in  Eq.  (9),  the  magnitude  of 

force  f-j  =  |f^.  I  has  the  following  form: 


4  =  24- 

cr 


cr 


r  S' 

a 


(11) 


o 


Figure  1.  Support  domain  of  particle  i 


3.  A  TWO-DIMENSIONAL  CASE  STUDY 


A  two-dimensional  case  is  considered  for  unit  cell 
nanocomposite  property  study.  The  case  is  similar  to  the 
case  in  reference  [15].  We  considered  that  the  epoxy 
resin  is  dispersed  well  by  5%  volume  of  nano  silicate 
particles.  The  particle  is  assumed  to  have  circular  shape 
to  simulate  the  configuration  of  the  nanoclay  flake.  The 
material  properties  of  silicate  are  given  as  below  [15]: 
density  2.2  g/cm^,  modulus  of  elasticity  73  GPa, 
Poisson’s  ratio  0.17,  tensile  strength  104  MPa.  The 
epoxy  properties  are  listed  as  below  [15]:  density  1.2 
g/cm^,  modulus  of  elasticity  2.3  GPa,  Poisson’s  ratio  0.4, 
tensile  strength  90  MPa.  In  the  epoxy  mixture  plate,  a 
square  plate  of  dimensions  40  x  40  x  0.05  pm^  is 
selected  as  a  typical  cell  model.  The  cell  model  is  fixed 
at  the  lower  edge  and  subjected  to  tensile  loading  of  60 
MPa  at  the  top  edge.  With  this  cell  model,  the  effects  of 
nanoparticle  size  on  the  mechanical  behavior  of  the 
material  were  studied. 

The  evenly  distributed  nanoparticle  case  was  studied 
first  with  the  regular  FEM  approach.  Five  cell  models 
are  considered  and  shown  in  Fig.  2.  Figure  2a  is  the  cell 
model  with  one  nanoparticle  of  10  pm  in  diameter.  The 
particle  size  is  obtained  by  calculating  the  requirement  of 
5%  volume  fraction  of  silicate  particles.  In  Figs.  2b-2e, 
the  particle  is  partitioned  to  4(2x2),  9(3x3),  36(6x6),  and 
144(12x12)  smaller  particles,  respectively,  with  the  same 
particle  volume  fraction.  The  corresponding  particle 
diameters  are  5,  3.4,  1.7,  and  0.84  pm.  Figure  2  also 


shows  von  Mises  stress  distribution  results  of  various 


cell  models  obtained  by  FEN 

1  simulations. 

# 

♦  ♦ 

♦  ♦ 

(a)  1  particle  (b)  4  particles 

♦  ♦  ♦ 


♦  ♦  ♦  .  .  *  *  . 

♦  ♦  ♦  I  I  I  I  I 

(c)  9  particles  (d)  36  particles 


(e)  144  particles 

Figure  2.  Cell  models  of  different  sizes  of  particles 
(FEM  approach) 

The  evenly  distributed  nanoparticle  case  was  then 
studied  with  the  meshfree  particle  method.  The  result  is 
compared  to  the  one  obtained  from  the  FEM  method.  In 
the  meshfree  particle  method,  the  epoxy  region  and 
nanoparticle  domain  are  discretized  into  small  particles. 
The  neighboring  particles  interact  with  each  other  in  a 
way  determined  by  the  potential  energy  between 
particles.  By  this  means,  no  particular  meshing 
technique  is  needed  for  particle  partitioning. 

Figure  3  shows  the  comparison  of  nanocomposite 
strength  and  stiffness  based  on  two  different  simulation 
methods.  The  variation  of  the  material  tensile  strength  is 
plotted  in  Fig.  3a  as  a  function  of  particle  dispersion.  As 
can  be  seen,  the  nanoclay-epoxy  composite  strength 
increases  as  the  nano  reinforcement  is  dispersed  using 
increasingly  fine  particles.  In  this  particular  example. 
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the  diameter  of  the  nanoelay  reaehes  840  nm  for  the 
partition  ease  144(12x12).  The  tensile  strength  of  the 
nanoeomposite  reaehes  its  peak  value  at  this  partition 
ease.  The  determination  of  material  strength  follows  the 
mle  suggested  by  Zhao  and  Hoa  [15]  in  the  study  of 
nanoelay-epoxy  eomposite  strength  for  a  similar  load 
ease.  Compared  to  the  FEM-based  result  (see  Fig.  3a), 
the  ehange  of  material  strength  is  more  apparent  in 
meshfree  partiele  method.  The  tensile  strength  obtained 
by  meshfree  partiele  method  eonverges  to  similar  value 
as  the  one  obtained  from  FEM  approach.  Consequently, 
the  precision  of  the  meshfree  particle  method  is 
acceptable  in  the  nanoeomposite  simulation. 

A  similar  tendency  can  also  be  seen  in  material 
tensile  modulus  (see  Fig.  3b),  although  the  change  of 
modulus  is  much  smaller  compared  to  the  change  in 
strength.  Furthermore,  the  tensile  modulus  does  not  vary 
significantly  as  the  function  of  particle  size,  which 
matches  the  reported  results  in  reference  [15]. 


In  Figs.  2  and  3,  the  nanoparticle  is  assumed  to  have 
same  size  and  to  be  distributed  evenly  in  the  epoxy  resin. 
In  practice,  it  is  natural  to  have  a  random  distribution  of 
particle  size  and  position  in  the  dispersion  process.  The 
effect  of  particle  randomization  is  studied  for  5%  volume 
fraction  of  nanoparticle  reinforcement  using  FEM 
approach.  Figure  4  shows  four  cell  models  and  stress 
plots  with  random  particle  distribution  obtained  from 
FEM  simulations.  The  corresponding  tensile  modulus 
and  strength  are  plotted  in  Fig.  5  and  are  compared  to  the 
results  for  evenly  distributed  cases.  As  can  be  seen,  the 
tensile  modulus  is  almost  identical  for  the  randomized 
cases  compared  to  the  evenly  distributed  ones,  while  the 
tensile  strength  will  depart  from  the  evenly  distributed 
counterpart.  The  difference  is  so  large  that  the  one 
particle  partition  case  will  have  the  best  tensile  strength 
after  randomization,  which  is  contradictory  to  the  evenly 
distributed  case  in  that  the  same  partition  has  the  worst 
strength  performance. 


(a)  tensile  strength 


(b)  tensile  modulus 

Figure  3.  Comparison  of  material  properties  based 
on  different  simulation  methods 
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(c)  9  particles  (d)  36  particles 

Figure  4.  Cell  models  for  nanoparticle 
randomization  study  (FEM  approach) 
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(a)  tensile  modulus 


(b)  tensile  strength 

Figure  5.  Comparison  of  tensile  modulus  and 
strength  for  randomized  and  even  nanoparticle 
distributions 


This  study  of  randomized  particle  distribution  shows 
the  importance  of  the  statistical  evaluation  for 
nanocomposite  properties  in  the  particle  dispersion 
process.  It  is  necessary  to  create  numerous  random 
samples  for  particle  partitions  under  given  particle  size. 
The  mechanical  properties  of  the  nanoclay-epoxy 
mixture  will  be  the  ensemble  averages  of  the  results  for 
all  samples.  The  statistical  tensile  strength  and  modulus 
have  been  studied  for  nanoclay-epoxy  composite  with  5 
volume  %  nanoparticle  concentration.  For  each  group  of 
particle  partition,  13  samples  were  created  to 
approximate  the  randomization  of  particle  size  and 
location.  Meshfree  particle  method  was  used  to  simulate 
the  material  property  for  each  sample. 


Particle  number 


(a)  tensile  modulus 


(b)  tensile  strength 

Figure  6.  Statistical  nanocomposite  property  for  5 
volume  %  nanoparticle  concentration 

Figure  6  shows  the  ensemble  averaged  value, 
maximum  value,  and  minimum  value  of  the  tensile 
modulus  and  tensile  strength  with  respect  to  various 
partition  groups.  As  can  be  seen,  the  change  of  tensile 
modulus  is  relatively  small  compared  to  the  variation  of 
tensile  strength.  For  partition  group  144  (the  diameter  of 
nanoparticle  is  about  840  nm),  the  minimum  tensile 
strength  is  approximately  two-thirds  of  the  maximum 
strength.  The  averaged  strength  is  closer  to  the 
minimum  value  than  to  the  maximum  one. 
Consequently,  statistical  analysis  is  cmcial  to  understand 
the  mechanical  properties  of  nanocomposites. 

In  the  final  numerical  example,  nanocomposite  with 
defects  was  considered  for  evaluating  tensile  strength. 
FEM  approach  is  used  to  simulate  the  nanocomposite 
property.  Figure  7  shows  the  stress  distribution  for  pure 
epoxy  and  nanoclay  embedded  epoxy.  The  defect  in  the 
composite  is  assumed  to  have  circular  shape  and  a 
diameter  close  to  1.0  micrometer.  The  circular  shaped 
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defect  is  used  to  mimic  the  bubble  created  during  the 
epoxy  manufacturing  process.  The  nanocomposite  has 
the  same  material  property  and  geometrical  configuration 
as  the  composite  shown  in  Fig.  2d.  The  epoxy  region  is 
evenly  distributed  with  36  nanoparticles,  while  the  9 
bubble  defects  are  also  evenly  distributed  in  the  epoxy. 
As  can  be  seen  in  Fig.  7,  the  maximum  stress  is  located 
near  the  edge  of  the  circular  shaped  defect.  The  tensile 
strength  for  nanoclay-epoxy  and  pure  epoxy  are 
compared  in  Fig.  7c.  The  tensile  strength  of 
nanocomposite  is  10  percent  higher  than  the  strength  of 
pure  epoxy.  As  a  result,  the  nanoclay  embedded  epoxy 
is  much  stronger  than  pure  epoxy.  Furthermore,  the 
strength  of  nanocomposite  is  relatively  insensitive  to  the 
existence  of  defect  compared  to  the  property  of  pure 
epoxy. 
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(c)  tensile  strength 

Figure  7.  von  Mises  stress  in  nanocomposite  and 
epoxy  with  defects  (FEM  approach) 

4.  NANOCLAY-REINFORCED  ARMOR  FOR 
BALLISTIC  PROTECTION 

MIOJfCIU  PMTIUTMe  PLAT! 

liM*  AfH 


I 
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(a)  nanoclay-epoxy  resin 

HtojKTU  PtrcnuriNO 


li 


(b)  pure  epoxy  resin 

Figure  8.  Penetration  process  for  composite  armor 
with  nanoclay-reinforced  and  pure  epoxy 

Developed  at  MKP  Inc.,  nanoclay-epoxy  resin  is 
used  as  the  gluing  polymer  in  the  ceramic  pellet 
faceplate  of  a  composite  armor.  It  has  been  found  that 
the  armor  with  nanoclay  insertion  maintained  the 
faceplate  integrity  after  two  rounds  of  shooting,  while 
the  faceplate  with  pure  epoxy  resin  damaged  after  the 
same  shooting  condition.  As  part  of  our  virtual 
prototyping  process,  a  numerical  simulation  was  utilized 
to  investigate  the  difference  in  protection  provided  by 
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armor  with  nanoclay  reinforced  epoxy  and  armor  with 
pure  epoxy  resin.  Efforts  have  been  made  to  reproduce 
similar  ballistic  protection  by  comparison  of  simulation 
to  actual  test  results.  In  this  part  of  research,  the 
improvement  of  mechanical  properties  for  epoxy  with 
nanoclay  was  represented  by  increases  of  modulus  and 
strength. 

The  simulation  was  based  on  LS-DYNA3D  and  the 
results  are  shown  in  Fig.  8.  It  can  be  seen  that  projectile 
(shown  in  red)  was  stopped  earlier  with  nanoclay- 
reinforced  matrix  material  (see  Fig.  8a).  For  armor  with 
pure  epoxy  as  the  bonding  media,  the  projectile  will 
penetrate  deeper  into  the  backplate  (see  Fig.  8b).  The 
improvement  of  armor  protection  can  be  explained  by 
the  fact  that  higher  strength  matrix  material  can  survive 
longer  and  contribute  more  to  the  defeat  of  projectiles. 

5.  CONCLUSIONS 

The  development  of  simulation  models  and  design 
methods  for  nanocomposite  will  enable  the  prediction, 
design,  and  prototyping  of  ballistic-protective  composite 
structures  in  a  broad  range  of  military  and  civilian 
applications  where  enhanced  structural  and  mechanical 
performance  is  required.  Those  that  potentially  benefit 
from  the  SBIR  program  include  DOD  organizations:  the 
Army,  Air  Force,  and  Navy.  Commercial  applications 
include  a  large  number  of  automotive  applications  to 
improve  the  strength,  stiffness,  weight,  surface  finish, 
flame  retardancy,  resistance  to  high  temperature  and 
corrosion,  among  many  desirable  characteristics. 
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